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<^ ! Abstract 

We give an elementary introduction to Hirota's direct method of constructing 
multisoliton solutions to integrable nonlinear evolution equations. We discuss in 
^ , detail how this works for equations in the Korteweg-de Vries class. We also show 

how Hirota's method can be used to search for new integrable evolution equations 
, and list the results that have been obtained before for the mKdV/sG and nlS classes. 

OO 
O 

. 

§ ; 1 Why the bihnear form ? 

In 1971 Hirota introduced a new direct method for constructing multisoliton solutions to 
^ integrable nonlinear evolution equations |I| . The idea was to make a transformation into 
■ new variables, so that in these new variables multisoliton solutions appear in a particularly 
simple form. The method turned out to be very effective and was quickly shown to give 
^ i A^-soliton solutions to the Korteweg-de Vries (KdV) modified Korteweg-de Vries 
(mKdV) 0, sine-Gordon (sG) ||^ and nonlinear Schrodinger (nlS) |^ equations. It is also 
useful in constructing their Backlund transformations |Q. Later it was observed that the 
new dependent variables (called "r-functions" ) have very good properties and this has 
become a starting point for further developments. 

In this talk our aim is to describe how multisoliton solutions can be constructed 
using Hirota's method. Multisoliton solutions can, of course, be derived by many other 
methods, e.g., by the inverse scattering transform (1ST) and various dressing methods. 
The advantage of Hirota's method over the others is that it is algebraic rather than 
analytic. The 1ST method is more powerful (it can handle general initial conditions) 
and at the same time more complicated. Accordingly, if one just wants to find soliton 
solutions, Hirota's method is the fastest in producing results. 
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2 From nonlinear to bilinear 



The (integrable) PDE's that appear in some particular (physical) problem are not usually 
in the best form for further analysis. For constructing soliton solutions the best form 
is Hirota's bilinear form (discussed below) and soliton solutions appear as polynomials 
of simple exponentials only in the corresponding new variables. The first problem we 
face is therefore to find the bilinearizing transformation. This is not algorithmic and can 
sometimes require the introduction of new dependent and sometimes even independent 
variables. 

Here we will discuss in detail only the KdV equation 

Uxxx + QuUx + Ut = 0. (1) 

Since we have not yet defined what bilinear is, let us first concentrate on transforming 
the equation into a form that is quadratic in the dependent variables. One guideline in 
searching for the transformation is that the leading derivative should go together with the 
nonlinear term, and, in particular, have the same number of derivatives. If we count a 
derivative with respect to x having degree 1, then to balance the first two terms of ([l|) u 
should have degree 2. Thus we introduce the transformation to a new dependent variable 
w (having degree 0) by 

u = dlw. (2) 
After this the KdV equation can be written 

'^xxxxx ~l~ ^^xx^xxx ~l~ '^xxt 0) (3) 

which can be integrated once with respect to x to give 

Wxxxx + ^wl^ + w^t = 0. (4) 

In principle this would introduce an integration constant (function of t), but since (^ 
defines w only up to w ^ w + x\(t), we can use this freedom to absorb it. 

Equations of this form can usually be bilinearized by introducing a new dependent 
variable whose natural degree (in the above sense) would be zero, e.g., logF or f/g. In 
this case the first one works, so let us define 

w = a log F (5) 

with a free parameter a. This results in an equation that is fourth degree in F, with the 
structure 

F^ X (something quadratic) + 3a(2 - a){2FF" - F'^)F'^ = 0. (6) 
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Thus we get a quadratic equation if we choose a = 2, and the result is 

FxxxxF — AFxxxFx + 3F^^ + FxtF — F^Ff = 0. (7) 

In addition to being quadratic in the dependent variables, an equation in the Hirota 
bilinear form must also satisfy a condition with respect to the derivatives: they should 
only appear in combinations that can be expressed using Hirota's D-operator, which is 
defined by: 

D:f-9 = {d..-dx,rf{xMx,)i^^^^^^. (8) 

Thus D operates on a product of two functions like the Leibnitz rule, except for a crucial 
sign difference. For example 

Dxf-g = fxg-fgx, 
F>xDt f-g = fgxt - fxgt - ftgx + fgxt- 

Using the D-operator we can write (|^) in the following condensed form 

{Dl + DxDt)F-F = 0. (9) 

To summarize: what we needed in order to obtain the bilinear form for (|l]) is a 
dependent variable transformation 

u = 2dllogF, (10) 

and we also had to integrate the equation once. 

For a further discussion of bilinearization, see e.g., P, |^. Unfortunately the process 
of bilinearization is far from being algorithmic. It is even difficult to find out beforehand 
how many new independent and/or dependent variables are needed for the bilinearization. 
In fact for some equations the natural form may not be bilinear but perhaps trilinear P] . 
Recently there have been some indications that singularity analysis can be used to reduce 
the guesswork: the number of dependent variables seems to be related to the number of 
singular manifolds 0. 

One important property of equations in Hirota's bilinear form is their gauge invariance. 
One can show 0] that for a quadratic expression homogeneous in the derivatives, i.e., of the 
form Y17=o {dl. f) g), the requirement of gauge invariance under / — > e''^ f, g — > e'^^g 
implies that the expression can be written in terms of Hirota derivatives. This gauge 
invariance can be taken as a starting point for further generalizations [§, p!0[] . 

Finally in this section we would like to list some useful properties of the bilinear 
derivative [Q. For P a polynomial, 

PiD)f-g = Pi-D)g-f, (11) 
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P{D)l-f = P{-d)f, 

dllogf 
5' log/ 



P{D)f-l = P{d)f, 
{Dlf-f)/{2n 

{Dtf.f)/{2f)-3{Dlf.ff/{2f). 



(12) 
(13) 
(14) 
(15) 



3 Constructing multisoliton solutions 

In this section we will construct soliton solutions for a class of equations. Properly speak- 
ing, the solutions should be called "solitary waves" until we can prove that they scatter 
elastically. However, since we are working towards integrable soliton equations, the soli- 
tary wave solutions will at some point become true solitons, so to make things simple we 
just call them solitons all the time. 



3.1 The vacuum, and the one-soUton solution 

Now that we have the KdV equation in the bilinear form, let us start constructing soliton 
solutions for it. In fact, it is equally easy to consider a whole class of bilinear equations 
of the form 

P{D.,,Dy,...)F-F = 0, (16) 

where P is a polynomial in the Hirota partial derivatives D. We may assume that P is 
even, because the odd terms cancel due to the antisymmetry of the D-operator. 

Let us start with the zero-soliton solution or the vacuum. We know that the KdV 
equation has a solution u = and now we want to find the corresponding F. From (|10[) 
we see that F = e20(t)x+/3(t) yig^g ^ that solves ([1|), and in view of the gauge freedom 
we can choose P = 1 as our vacuum solution. It solves (|16D provided that 

P(0,0,...) =0. (17) 

This is then the first condition that we have to impose on the polynomial P in ([16|) . 

The multisoliton solutions are obtained by finite perturbation expansions around the 
vacuum P = 1: 

P = l + e/i + eV2 + e'/3 + --- (18) 

Here e is a formal expansion parameter. For the one-soliton solution (ISS) only one term 
is needed. If we substitute 

F = l + eh (19) 

into ( [TB| ) we obtain 

P(D,,...){l-l + el-/i + e/i-l + eVi-/i} = 0. 
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The term of order e° vanishes because of ([T7|) . For the terms of order we use property 
( [T^ ) so that, since now P is even, we get 

P{d,,dy,...)f,=0. (20) 

The sohton solutions correspond to the exponential solutions of For a ISS we take 
an fi with just one exponential 

fi = e"^, r] = px + qy + . . . + const, (21) 

and then ( pO]) becomes the dispersion relation on the parameters p,q, . . . 

P{p,q,...) = 0. (22) 

Finally, the order term vanishes because 

P(D)e''-e'' = e^''P{p-p) = 0, 

by (0). 

In summary, the ISS is given by (|l^,^) where the parameters are constrained by (p2|). 
For KdV the dispersion relation is = p. 

3.2 The two-soliton solution 

The 2SS is built from two ISS's, and one important principle is that for integrable systems 
one must be able to combine any pair of ISS's built on top of the same vacuum. Thus if 
we have two ISS's, Fi = 1 + e^^ and F2 = l + e"^'^, we should be able to combine them into 
a form F = 1 + /i + /2, where /i = e''^ + e''^ Gauge invariance suggest that we should 
try the combination 

F = 1 + e"^ + e"-' + ^126"^+^"^ (23) 
where there is just one arbitrary constant A12. Substituting this into ([TBI) yields 



P{D){ 1 


1 


+ 


1 ■ e'^i 


+ 


1 ■ e''2 


+ 




1 


+ 


gr,! . gr,i 


+ 


gm . g»?2 


+ 




1 


+ 


gr?2 . gr?i 


+ 


g'72 . gr?2 


+ 




1 


+ 


Al2e''i+''2 . gm 


+ 


Al2e'^^+''^^ ■ e"^^ 


+ 



In this equation all non- underlined terms vanish due to ( p!7| , p2[ ). Since P is even, the 
underlined terms combine as 2Ai2P{pi + P2) + 2P(pi — P2) = 0, from which A12 can be 
solved as 

A,2 = -^^^^. (24) 
P{pi + P2) 

The important thing about this result is that we were able to construct a two-soliton 
solution for a huge class of equations, namely all those whose bilinear form is of type (pIBI). 
In particular this includes many non-integrable systems. 
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3.3 Multi-soliton solutions 



The above shows that for the KdV class ([Tq) the existence of 2SS is not strongly related 
to integrability, but it turns out that the existence on 3SS is very restrictive. 

A 3SS should start with /i = e^^ + e"^^ + e^^ and, if the above is any guide, contain 
terms up to f^. If we now use the requirement that the solution should reduce to a 2SS 
when the third soliton goes to infinity (which corresponds to 77^ — > ±00) then one finds 
that F must have the form 

F = 1 + e''! + e''^ + e''^ 

+^126"^+"^ + Ai3e''i+''^ + ^236''^+''^ + ^12^13^236"^+"^+''^ (25) 

Note in particular that this expression contains no additional freedom. The parameters 



Pi are only required to satisfy the dispersion relation (|2^) and the phase factors A were 



already determined (plj). This extends to NSS |n 



N 



F= ^ exp 'Pihj)^il^j + , (26) 



Mi=0,i \l<j<7<Af i=l 

l<i<N ^ 



where {Aij = e'^'-*'-'^). Thus the ansatz for a NSS is completely fixed and the requirement 
that it be a solution of (|16D implies conditions on the equation itself. Only for integrable 
equations can we combine solitons in this simple way. More precisely, let us make the 

Definition: A set of equations written in the Hirota bilinear form is Hirota inte- 
grable, if one can combine any number N of one-soliton solutions into an NSS, and the 
combination is a finite polynomial in the 's involved. 

In all cases known so far, Hirota integrability has turned out to be equivalent to more 
conventional definitions of integrability. 



4 Searching for integrable evolution equations. 

Since the existence of a 3SS is very restrictive, one can use it as a method for searching 
for new integrable equations. All search methods depend on some initial assumptions 
about the structure. In this case we assume that the nonlinear PDE can be put into a 
bilinear form of type (|TB|), but no assumptions are made for example on the number of 
independent variables. 
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4.1 KdV 



If one now substitutes (|2|) into (|16D one obtains the condition 



(27) 



where the symbol = means that the equation is required to hold only when the param- 
eters Pi satisfy the dispersion relation P{pi) = 0. 

In order to find possible solutions of (^7|) we made a computer assisted study [|12| and 
the result was that the only genuinely nonlinear equations that solved ( P?]) were 



(D^ - AD,Dt + 3Dj)F ■ F 
{DID, + aDl + D,Dy)F ■ F 
{Dl - D^Dl + aDl + bD^Dt + cDl)F ■ F 
{Dl + hDlDt - hDl + D,Dy)F ■ F 



0, 
0, 
0, 
0. 



(28) 
(29) 
(30) 
(31) 



and their reductions. These equations also have 4SS and they all pass the Painleve 
test |]14|. Among them we recognize the Kadomtsev-Petviashvili (containing KdV and 
Boussinesq) (pH), Hirota-Satsuma-Ito (E^) and Sawada-Kotera-Ramani (|3TD equations; 



they appear in the Jimbo-Miwa classification [0. The only new equation is (|30D. It is 
somewhat mysterious. It has not been identified within the Jimbo-Miwa classification 
because it has no nontrivial scaling invariances, furthermore we do not know its Lax pair 
or Backlund transformation. But it does have at least 4SS, and it passes the Painleve 
test PI. 



4.2 mKdV and sG 

As was mentioned before, Hirota's bilinear method has been applied to many other equa- 
tions beside KdV. Here we would like to mention briefly some of them. 

For example the modified Korteweg-de Vries (mKdV) and sine-Gordon (sG) equations 
have a bilinear form of the type 

BmG-F = ,, (32) 
A(Di)(F-F + G-G) = li, 

where A is even and B either odd (mKdV) or even (sG). For mKdV we have B = 
Dl + Dt, A = Dl, and for SG, B = D^Dy — 1, A = D^Dy. This class of equations also has 
2SS for any choice of A and B. If B is odd one can make a rotation F = f+g, G = i{f~g) 
after which the pair ( ^2]) becomes B g ■ f = 0, A g ■ f = 0. 
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In principle the pair (|32[) can have two different kinds of sohtons, 

= 1 + e^^, G = 0, with dispersion relation A{p) = 0, 
= 1, G = e''^, with dispersion relation B{p) = 0. 

In mKdV and SG the A polynomial is too trivial to make the first kind of soliton inter- 
esting. In [|T5], 1^ we searched for polynomials A and B for which any set of three solitons 



could be combined for a 3SS. The final result contains 5 equation, of mKdV type, three of 
them have a nonlinear B polynomial but a factorizable A part (and hence only one kind 
of soliton with B acting as the dispersion relation), 

(34) 
(35) 
(36) 



{aDl + bDl + DlDt + Dy)G- 


F = 


0, 


DIG- 


F = 


0, 


{ {aDl + hDl + Dy)G- 


F = 


0, 


[ D.DtG- 


F = 


0, 


f {D,DyDt + aD, + bDt)G 


■F = 


0, 


\ D^DtG 


■F = 


0. 



For a discussion of the nonlinear versions of the last two equations, see |T6l . 

In two cases both A and B are nonlinear enough to support solitons, note that the B 
polynomials are the same and that the A parts have already appeared in the KdV list: 



{Dl + Dy)G-F = 0, 

{DlDt + aDl + DtDy) G ■ F = 0, 

{Dl + Dy)G-F = 0, 
Two equations of sine-Gordon type were also found: 



{Dl + 5DlDy-5Dl + DtD.,)G-F = 0. 



(37) 
(38) 



{D,Dt + b)G-F = 0, 

{DlDt + 3bDl + DtDy){F ■ F + G ■ G) = 0, 

{aDlDt + DtDy + b) G ■ F = 0, 

D^Dt{F ■ F + G ■ G) = 0. 



(39) 
(40) 



4.3 nlS 



A similar search was performed [0, |^, [T^ on equations of nonlinear Schrodinger (nlS) 
type, 

' B{D,)G.F = 0, 
A{D,)F-F = ^ ' 
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f (Dl + iDy + c)G- 


■F = 


0, 


\ {a{Di - 3Dl) + D^Dt) F ■ 


■F = 


|G|2 


{iaDl + 3cL)2 + _ 2dDt) + g) G ■ 


■F = 


0, 


aD'lDt + aDl + [b + 3c^)D^Dt + dD^) F • 


■F = 


|G|2 


J {iaDl + 3D^Dy - 2iDt + c) G ■ 


■F = 


0, 


\ {a{a''Dl - + AaD^Dt) + hDl) F ■ 


■F = 





where F is real and G complex. Again two kinds of solitons exist, 

F = 1 + e''-^, G = 0, with dispersion relation A^p) = 0, ^^^^ 

F = 1 + Ke^'^^'^B ^ G = e''^, with dispersion relation B{p) = 0. 

In this case the existence of a 2SS is not automatic because the ISS already involves terms 
of order and a 2SS therefore e^, whereas in the previous cases contributions were 
sufficient for 2SS. 

Three equations were found in this search: 

(43) 
(44) 
(45) 

Perhaps the most interesting new equation above is the combination in ( ^Sf ) of the 
two most important (2 + l)-dimensional equations, Davey-Stewartson and Kadomtsev- 
Petviashvili equations, see [jl8|, page 315. 

The above lists contains many equations of which nothing else is known other than that 
they have 3SS and 4SS. This alone suggests that they are good candidates for integrable 
2 + 1 dimensional equations, but connections to other definitions of integrability (like Lax 
pairs) are still open. 
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